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We propose a practical scheme to generate a pure valley current in monolayer transition metal 
dichalcogenides by one-photon absorption of linearly polarized light. We show that the pure valley 
current can be detected by either photoluminescence measurements or the ultrafast pump-probe 
technique. Our method, together with the previously demonstrated generation of valley polarization, 
opens up the exciting possibility of ultrafast optical-only manipulation of the valley index. The tilted 
field effect on the valley current in experiment is also discussed. 
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I. INTRODUCTION 


Recent years have seen a surge of interest in the manip¬ 
ulation of the valley index of Bloch electronslargely 
driven by its potential applications in electronics and op¬ 
toelectronics.^ The valley index enumerates degenerate 
energy extrema in momentum space. Such degeneracies 
are often present in 2D materials with a honeycomb-like 
structure, such as graphene, boron nitride, and transi¬ 
tion metal dichalcogenides (TMD). In these materials, 
weak intervalley scattering renders the valley index an 
effective degree of freedom that can be utilized in novel 
devices. This realization of valley-based electronics is 
called valleytronics, which depends crucially on the dy¬ 
namical control of two quantities: valley polarization 
and valley current. The optical gener ation of valley po¬ 
larization by circularly polarized lightP^E] shows some 
promise as a realization of valleytronics and has been ex¬ 
perimentally demonstrated in monolayer On 

the other hand, so far only valley-polarized electric cur¬ 
rent has been reported.^i^*l^In analogy to spintronics, it 
would be desirable to generate a pure valley current, in 
which there is no net motion of charge; carriers in oppo¬ 
site valley move in opposite direction. Such a pure valley 
current would rule out any charge-related effedP®^ and 
generate minimal Joule heating, similar to a pure spin 
currenlP^. 

In this work, we propose a new approach to the gen¬ 
eration and detection of a pure valley current by optical 
means. Based on both symmetry analysis and an effec¬ 
tive k p Hamiltonian, we show that a pure valley current 
can be generated by linearly polarized light in monolayer 
TMDs. The generating mechanism parallels that for spin 
current.^ However, the role of spin-orbit coupling is re¬ 
placed by the trigonal warping in the band structure, 
which is entirely a lattice effect. Due to the unique spin- 
valley coupling in this system,li^the generated valley cur¬ 
rent is accompanied by a spin current. We also present 
a theory for valley diffusion that takes into account the 
spin-valley coupling, and show that the pure valley cur¬ 
rent can be detected by either photoluminescence mea¬ 
surements or the ultrafast pump-probe technique. Our 


method, together with the previously demonstrated gen¬ 
eration of valley polarization,lHHi^ opens up the exciting 
possibility of ultrafast optical-only manipulation of the 
valley index. 

The paper is organized as follows. In Sec. |TTj we 
present the optical generation of valley current, where the 
symmetry analysis is given in Sec. |H A| and the numerical 
result is shown in Sec. 11 Detection of the generated 
valley current is considered in Sec. m where the ultra¬ 
fast pump-probe and photoluminescence measurements 
are proposed in Sec. HI A| and Sec. HIB respectively. 
Finally, discussion and conclusion are made in Sec. |IV| 


II. OPTICAL GENERATION OF VALLEY 
CURRENT 

A. Symmetry analysis 

Figure shows the schematic setup. A linearly po¬ 
larized light E{uj) = Eq{cos9x sin 6y) at normal in¬ 
cidence is considered, where Eq and 9 refer to the am¬ 
plitude and polarization angle of the electric field, re¬ 
spectively. We choose x to be along the zigzag direction 
and y the armchair direction. In monolayer TMD, each 
transition metal cation is trigonal-prismatically coordi¬ 
nated by six nearest neighbor chalcogen anions, explicitly 
breaking the inversion symmetry. The relevant symme¬ 
try operations of the system include three-fold rotation 
Ca around the z axis, mirror reflections Mx{x —^ —x) 
and Mz{z —>■ —z), and time-reversal. 

The generation of a dc current by imposing an op¬ 
tical field—namely, the photogalvanic effect (PGE)—is 
a second-order nonlinear phenomenon characteristic of 
non-centrosymmetric materials. Under a monochromatic 
light E(t) = E{uj)e~‘^‘^* c.c., the photocurrent has the 

standard form 

3a=^XaM{k,uj,-uj)Ep{uj)E*^{uj) , 
k 

where the fe-resolved second-order susceptibility tensor 
Xai 3 -y{k,uj,—uj) is given from the perturbation theory 
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FIG. 1. (color online). The crystal structure of monolayer 
TMD, where M (red) is the transition metal atom and X 
(blue) is the chalcogen atom, (a) Side and (b) top view of the 
lattice structure. The thick arrow depicts a normally incident 
linearly polarized light with the polarization angle 9. The yz 
plane is dehned as one of the mirror planes (shaded region). 


bjEmn 


Xa0j{k,uj,-uj) = 2^ f dti ( dt2e 

n J-oo J-oo 

Xc,/3-f(k,ti,t2) = TY{po{k)[[Va,Vi3{ti)],V^{t2)]). 

( 2 ) 

Po{k) is the initial equilibrium density matrix operator 
and Va,i 3 ,~fit) are the velocity operator in the Heisenberg 
picture at time t. Tr denotes the trace and S is the area. 
Since we are dealing with a strictly 2D system, the in¬ 
dices a, /3 ,7 can be either x or y. According to Ref. 
Xa/3-y{k,uj, —uj) is further decomposed into three terms: 
two constant terms and one linear-in-time term, the lat¬ 
ter of which is cut off by relaxation time r based on the 
relaxation time approximation. In this work, we restrict 
ourselves to the high quality samples with t ^ h/AE, 
in which the response is dominated by the linear-in-time 
term, and the other two terms can be neglected.!^ AE is 
the optical transtion gap. As a result, the susceptibility 
reduces to 


Uj) — ^ (^/3,fc)nm('C7,/g)mr 

n,m 

X {,’^a^k')nn\d^nm^k 

X 5(J%UJ Sjyik ^nk) ■ 


(3) 


Here {va,k)mn = {nT-, k\va\n, k) is the velocity matrix ele¬ 
ment in the Bloch basis |n, fc), is the band dispersion, 

d^nm.k — fnk fnik with fnk “ [1 “b 6Xp[/3(ey2^/5, M)]] 

being the Fermi-Dirac distribution, p is the chemical po¬ 
tential and /3 = l/ksT. Note that the expression for ja, 
when transformed into the real spac e repr esentation, is 
also recognized as the “shift current” IHIIZI 

For a system with time-reversal symmetry (TRS), 
PGE vanishes under linearly polarized light. The rea¬ 
son is that linear polarization picks out the real part of 



FIG. 2. (color online). The first Brillouin zone (BZ) of mono- 
layer MX 2 . Pi, 2 , 3,4 are symmetry-related points. 


X in Eq. ([^, which satisfies 

xlp^ik) = -xip-yi-k) , xip^ik) = -xlp-yi-k) (4) 

due to the TRS. For simplicity, we have omitted the ar¬ 
guments UJ and —UJ in x- Summing over each pair of x’s 
in Eq. @ then yields zero charge current. Equation 0 
suggests that it is possible to generate a pure spin cur¬ 
rent. However, without breaking the spin degeneracy, 
one has xl^p^ik) — Xap-y{k)', consequently, the total spin 
current still vanishes. In Ref. mi Bhat et al. showed that 
introducing the spin-orbit coupling can break the spin 
degeneracy, giving rise to a pure spin current. 

We now show that a similar effect can generate a pure 
valley current, i.e., -I- = 0, and yb 0. In mono- 

layer TMD, the two valleys, located at the K and K' 
points of the hexagonal Brillouin zone, are related by the 
TRS. As such, the valley-resolved susceptibility tensor 
satisfies 

Xap-yi(l) = -X%-y{-(l) , X%-y{.(l) = -XaPji-g) , (5) 

where q = k — K{K') defines a small momentum away 
from the valley center K {K'). Again, the charge current 
vanishes. However, the K and K' points have Ca rota¬ 
tional symmetry. This allows Xap-yil) 7 ^ Xapji^)^ lead¬ 
ing to a pure valley current oc '£,qJ2r,^vXap-yi9)’^^^'^^ 
= ± 1 for 77 = K, K'. For a system symmetric un¬ 
der Mz, C3 is the only possible rotation symmetry that 
can break the valley “degeneracy” and induce the valley 
current; all other rotational symmetries yield zero valley 
current, even though they given rise to an anisotropic 
band structure. 

To substantiate the preceeding argument, we carry out 
a detailed group theory analysis. The symmetry group 
of the K point is Csh', X transforms as a direct product 
E'®E'®E', which contains two copies of the identity rep¬ 
resentation. This indicates that there are two indepen¬ 
dent components of y: X^yy = “XLy = -Xlyx = -X^xx 
(denoted by xT) and xlxx = -Xxyy = -X^xy = -Xyyx 
(denoted by x2)- However, due to the TRS, we find that 
the combined symmetry TM^ requires the contribution 
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from Xe to vanish. To see this, let us consider two k- 
points Pi and P 2 in the K valley. These are related by 
Qy —% We also introduce two intermediate 

points P 3 and Pa. Pi and P 3 are related by the mirror 
symmetry and satisfy 

Xei(lx,qv)=xf(.-qx,qy)- ( 6 ) 

Meanwhile, P 3 and P 4 are related by the TRS, so that 

Xeiqx,qy) = -Xe ' {-qx,-qy) ■ (7) 

Finally, P 2 and P 4 are equivalent up to a reciprocal lattice 
vector. Therefore, upon summing over fe-states in one 
valley, the contribution from x'l vanishes. The only non¬ 
vanishing contribution to the valley current is from 
with an angular dependence 



Photon energy - band gap (eV) 


f oc xlxx cos(26» -P ip) , (8) 

where ip is the detection angle. Experimentally, by fixing 
(p, a TT-period oscillation of the signal is expected. 


B. Numerical results 


A unique property of monolayer TMD is the strong 
spin-valley coupling, which refers to the opposite spin 
splitting at the valence band edge in opposite valleys [see 
Fig.[^a)].li^We can see immediately that a pure spin cur¬ 
rent will accompany the generated valley current. Note 
that this spin current would vanish if the energy bands at 
the K and K' points are isotropic, even in the presence 
of spin-orbit coupling. 

Another important parameter in Eq. is the relax¬ 
ation time T. Due to its multivalleyed band structure, 
there are several scattering channels in monolayer TMD. 
They give rise to four relaxation times '^li^.^alinter ’'^Tiich 
refer to intra- and intervalley scatte ring by electrons (e) 
and holes {h). These lifetimes satisfjBSHani 


T T T' 

intra inter 


h 

inter 


(9) 


Given that l/r = 1 /Tintra + 1 !Winter , the optically gener¬ 
ated valley current is predominantly determined by the 
intra-valley scattering time Tintra- Electrons and holes 
exhibit almost the same intra-valley scattering time when 
neglecting the weak intrinsic electron-hole asymmetry in 
the system. The last inequality in Eq. ® comes from 
the aforementioned spin-valley coupling Although 
not essential in the generation of the valley current, it 
is important to the detection process as discussed below. 
In the absence of spin-dependent scattering, the upper 
limit of is set by the Bir-Ar onov- Pikus mechanism, 
which could be as large as ~ Ins.^^^ISIl 

To calculate the valley current in monolayer TMD, we 
employ a low-energy effective k ■ p HamiltoniaiP^ around 
valley K{K'), which includes both the C 3 symmetry- 
allowed trigonal warping and fc-cubed corrections. We 


FIG. 3. (color online). Linear valley and spin photogalvanic 
effect (PGE) of monolayer M 0 S 2 . (a) Schematics of band dis¬ 
persion around valley K{K') along the axis. Red (green) 
curves denote states with spin up (down). Thick (thin) solid 
blue arrows depict strong (weak) optical transition rates in 
each valley. Fermi energy Ep = 0. (b) Angular dependence 
of valley current on the polarization angle 9 along the x axis 
(zigzag direction). Red solid and blue dashed curves label 
photon energy Tiuj = 1.68eV and l.SOeV. (c) Valley and spin 
current as functions of photon energy (minus by band gap 
Eg), r = 55 fs, A = 3.01 x 10'* V/m, T = 5 K. Band param¬ 
eters are adopted from Ref. |3^and band gap Eg = 1.585eV. 


also take into account the large spin splitting in the va¬ 
lence bands; the small spin splitting in the conduction 
bands is ignored. 


With realistic parameters,^ our numerical results are 
shown in Fig. In Fig. [^a), the band dispersion of 
the effective model is plotted for ky = 0 , which clearly 
shows the large spin splitting in the valence bands. Due 
to the symmetry, optical transitions excite states with 
different \qx\ in each valley. This results in different op¬ 
tical transition rates (indicated by the thickness of the 
arrow) and different velocity in the Qx direction, both of 
which contribute to generating the valley photocurrent. 
Figure l^b) shows the angular dependence of the valley 
current on the polarization angle 6 by fixing (p = Q. The 
valley current j'" = jk — jx' has an order of 10“® A/cm, 
comparable to the magnitude of photocurrent observed in 
GaAs quantum wells.l^ Figure [^c) displays both valley 
and spin current as functions of photon energy. We note 
that as soon as the lower spin-split valence band becomes 
active, the spin current displays a downward trend. This 
allows us to manipulate the generation of valley and spin 
current either collectively or separately. 
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III. DETECTION OF VALLEY CURRENT 

Next we discuss the detection of the pure valley cur¬ 
rent. Our idea utilizes the fact that the valley carriers in 
monolayer TMD are described by a pair of massive Dirac 
fermions with opposite mass,^ therefore each valley ex¬ 
hibits opposite time-reversal symmetry breakin g ef fects 
such as circular dichroisrrfi^ and Faraday rotationl^Note 
that there is a possible complication due to the lar ge e x- 
citon binding energy observed in monolayer TMDpSEo] 
which makes the generation of free carriers difficult. To 
remove the exciton effect, we may heavily dope the sam- 
pl^^ or apply a large in-plane electric field.^^ Under 
these circumstances, we propose two possible detection 
schemes. 


A. Ultrafast pump-probe measurement 

In one scheme, one can observe the second-harmonic 
generation (SHG) of the valley current using the ultra¬ 
fast pump-probe techniqu^^as shown in Fig. [^a). This 
is similar to the proposed detection method of spin cur- 
rent.l^ First, the pump light (with frequency /i > Eg) 
generates a pure valley current. Then the probe light 
(with frequency /2 < Egj2) creates a population im¬ 
balance between k and —k states, which leads to a 
net Faraday rotation. This will induce a polarized field 
P(2/2) oc {z X E\E\), and emit a second-harmonic signal 
(orthogonal to the probe light). Since the energy of the 
SHG is still below the band gap Eg, the Faraday rota¬ 
tion is related to the virtual interband transition, which 
distinguishes it from other optical effects of the pump 
light. 


B. Photoluminescence measurement 

An alternative proposal, unique to monolayer TMD, 
is to investigate the photoluminescence (PL) helicity as 
shown in Fig. ib). Suppose a linearly polarized light 
illuminates the central region, generating a steady valley 
current. As the valley carriers move outside of the central 
region, they will start the diffusion process described by 

DV^i5/i(a;) - 5^i{x)/nnter = 0 , (10) 

where D = v\t j2 is the diffusion constant—derived from 
the Fermi velocity vp and the momentum relaxation time 
T—and Sfj,{x) = ijl^{x) — (x) is the chemical poten¬ 

tial difference between the two valleys. This equation 
describes both electrons and holes. Gonsider the right re¬ 
gion of the sample. For a valley current with initial veloc¬ 
ity j'’ = ^{<^xx/2e)dx6n{x), where -l-(—) corresponds to 
holes (electrons) and axx is the total longitudinal (Drude) 
conductivity, we obtain 

(5/i(x) = (T^CJ iinter/^xx) X/£inter) , (H) 


(a) 



2f2 



FIG. 4. (color online). Detection of pure valley current in 
monolayer M 0 S 2 . (a) The second-harmonic generation from 
the pump-probe experiment. Red (green) line denotes pump 
(probe) light, and purple line refers to the second-harmonic 
generation. /i(/ 2 ) is the frequency of the pump (probe) light, 
(b) Photoluminescence (PL) behavior over the electron-hole 
recombination length d. Red (blue) ball labels carriers from 
valley K(K') and ’-’(’-!-’) denotes electrons (holes). 


where linter = VDTinter IS the valley diffusion length. In 
monolayer TMD, we have <C Team ^ , where 

Tcom is the electron-hole recombination time.^ There¬ 
fore, after a diffusion length d ~ V Dreom : holes will 
have a local chemical potential difference Sfi^ix) < 0; 
meanwhile, electrons become almost equally populated 
in the two valleys and 5^e{x) ~ 0. Following the valley¬ 
contrasting circular dichroism,E2HIII this leads to a net 
a~*' PL hecility. A similar argument can be applied to 
the left region of the sample, where a net cr“ PL is ex¬ 
pected. Experimental results founeP^ that ~ I 

ns, which leads to €^„(gr ~ 1 /rm; they further deter¬ 
mined that Team = 50 ps, resulting in a PL helicity 
(p^ — )/{e^ + ) = 0.8. Note that this phe¬ 

nomenon is intimately related to the spin-valley coupled 
bands, and hence absent in other multi-valleyed systems 
such as staggered monolayer or biased bilayer graphene. 


IV. DISCUSSION AND CONCLUSION 

So far we have considered the normal incidence case. 
For oblique incidence, the results are summarized in Ta¬ 
ble m where the general form of the valley-resolved cur¬ 
rent jxijy) along the longitudinal (transverse) direction 
is given, rj refers to valleys K and K', and the xz plane 
is the incident plane. For a linearly polarized light, val¬ 
ley current is induced in both longitudinal and transverse 
directions; in contrast, only the longitudinal valley cur¬ 
rent is generated under circularly polarized light. In both 
cases, charge current vanishes since Xxxx + Xxxx = 0 due 
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TABLE I. Valley-resolved longitudinal (transverse) current j2{jy) under oblique incidence for linear and circular photogalvanic 
and photon drag effect, with rj = K ox K'. 9 and 4> are polarization and incident angle, respectively. T:^xxx, T^xyy and T^xxy 
are independent components of a rank-4 tensor. Eq is the amplitude of the electric field. 



Longitudinal {j'2) 

Transverse (j()) 

Linear photogalvanic effect 

X^xxicos'^ G cos'^ (j) - sin^ 9) El 

-XxxxCos(l>sm 29El 

Circular photogalvanic effect 

-Xxxx sin^ (j>El 

0 

Linear photon drag effect 

Qx {TSx^xx cos^ 9 cos^ 4> + TS^yy sin^ 9)El 

% {T^xxx - T^xyy) COS (f) siu 29EI 

Circular photon drag effect 

Qx{TSxxx cos^ (j) + TS^yy - 2iTSxxy cos (i))El 

0 


to the mirror symmetry under leading to a pure val¬ 
ley current. Apart from the PGE, there exists another 
photocurrent generating mechanism under the oblique 
incidence, namely the photon drag effect (PDE).I^ In 
this case, photons transfer both momentum and angu¬ 
lar momentum to carriers, and the current is described 
by ja = where T is a rank-4 tensor and 

Q is the photon wavevector. Similar to the PGE, PDE 
contributes to the valley current under both linearly and 
circularly polarized light. However, ^ 0, 

T^xxyyET^^yy ^ 0, indicating that the net charge current 
does not necessarily vanish in the system. To distin¬ 
guish these two mechanisms in experiments, one notices 
the fact that the response from PDE (PGE) is an odd 
(even) function under the reversal of incident direction 
Q —>■ —Q, by which the dominant mechanism can be 
identified. 

The proposal of pure valley current generation can be 
generalized to other systems with appropriate symme¬ 
tries, however, the magnitude of the effect and the detec¬ 
tion scheme may vary among different systems. In our 
view, monolayer TMD has its advantage that the spe¬ 


cific band and symmetries allow peculiar detection and 
observable signal. 

To conclude, we have demonstrated that a linearly- 
polarized light can induce a pure valley current in mono- 
layer TMD. This mechanism originates from the C^- 
symmetry rather than spin-orbit coupling. Furthermore, 
we propose two realistic optical approaches to detect the 
pure valley current. The effect of oblique incidence is also 
discussed. 

Note added .—Upon the completion of this work, we 
become aware of two recent papers, Ref. |46] and 133 which 
also studied the nonlinear valley effect. 
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